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ABSTRACT 

We present the results of a Chandra observation of the galaxy cluster Abell 689 
(z=0.279). Abell 689 is one of the most luminous clusters detected in the ROSAT 
All Sky Survey (RASS), but was flagged as possibly including significant point source 
contamination. The small PSF of the Chandra telescope allows us to confirm this and 
separate the point source from the extended cluster X-ray emission. For the cluster we 
determine a bolometric luminosity of Lboi = (3.3±0.3) x 10 44 erg s _1 and a temperature 
of kT=5.1^^3 keV when including a physically motivated background model. We 
compare our measured luminosity for A689 to that quoted in the Rosat All Sky Survey 
(RASS) and find Lo.i-2.4.kcV = 2.8 xlO 44 erg s _1 , a value ~10 times lower than the 
ROSAT measurement. Our analysis of the point source shows evidence for significant 
pileup, with a pile-up fraction of ~60%. SDSS spectra and HST images lead us to 
the conclusion that the point source within Abell 689 is a BL Lac object. Using radio 
and optical observations from the VLA and HST archives, we determine a ro =0.50, 
«ox=0.77 and a rx =0.58 for the BL Lac, which would classify it as being of 'High- 
energy peak BL Lac' (HBL) type. Spectra extracted of A689 show a hard X-ray 
excess at energies above 6 keV that we interpret as inverse Compton emission from 
aged electrons that may have been transported into the cluster from the BL Lac. 
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1 INTRODUCTION 

Studies of clusters of galaxies, including measurements of 
their number density and growth from the highest density 
perturbations in the early Universe, offer i nsight into the un - 
derlying cosmology (|Vikhlinin et alj2003l ; [Allen et alj2004h . 
However, in order to use clusters as a cosmological prob e 
three essential tools are required (|Del Popolo et al.ll20ich : 
(a) an efficient method to find clusters over a wide redshift 
range, (b) an observational method of determining the clus- 
ter mass, and (c) a method to compute the selection function 
or the survey volume in which clusters are found. These re- 
quirements are met by large surveys with well understood 
selection criteria. Arguably the most effective method of 
building large, well defined cluster samples has been via X- 
ray selection. The high X-ray luminosities of clusters make 
it relatively easy to detect and study clusters out to high 
redshifts. 

Many cluster samples have been constructed based 



E-mail: P.Giles@bristol.ac.uk 



upon large X-ray surve ys such as the Eins tein Medium Sen- 
sitivity Survey (EMSS; iGioia et alJll99Ch and the ROSAT 
All Sky Survey (RASS: FVoges! Il992ft . However due to the 
relatively poor angular resolution of these X-ray observato- 
ries, observations of clusters were susceptible to point source 
contaminatio n. Indeed, within the ROSAT Brightest Cluster 
Sample (BCS; Ebeling et aljl998h and its low-flux extension 
(eBCS: lEbeling et alj|200cfc 9 out of 201 clusters and 8 out 
of 99 clusters respectively were flagged as probably having a 
significant fraction of the quoted flux from embedded point 
sources. Being able to resolve these point sources is of crucial 
importance for the reliable estimation of cluster properties, 
and indeed the nature of the point source contamination is 
of independent interest. 

The study of galaxy clusters has been transformed with 
the launch of powerful X-ray telescopes such as Chandra and 
XMM Newton, which have allowed the study of the X-ray 
emitting intracluster medium (ICM) with unprecedented de- 
tail and accuracy. With the launch of this new generation 
of X-ray telescopes, we are able to uncover interesting fea- 
tures in the morphologies of individual clusters. In particu- 
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lar, Chandra's high angular resolution provides the means 
to examine individual cluster feat ures with gr eat detail. 

Abell 689 (her eafter A689; lAbelll Il958h . at z=0.279 
l|Collins et all 1 19951 ). was detected in the RASS in an ac- 
cumulated exposure time of 317s. It is included in the BCS, 
with a measured X-ray luminosity of 3.0xl0 45 erg s _1 in 
the 0.1-2.4 keV band. This luminosity is the third high- 
est in the BCS, and thus A689 meets the selection crite- 
ria for various hi ghly X-ray luminous cluster samples (e.g. 
iDahle et al]|2002h . However this cluster was noted as having 
possible point source contamination, and for this reason has 
often been rejected from many flux limited samples. In this 
study we present results of Chandra observations designed to 
separate any point sources and determine uncontaminated 
cluster properties. 

The outline of this paper is as follows. In §[2] we discuss 
the observation and data analysis. Results of the X-ray clus- 
ter analysis is presented in §[3] In §[4]we present our analysis 
of the X-ray point source through X-ray, optical and radio 
observations. We interpret our results in §[5]and the conclu- 
sions are presented in § [6] Throughout this paper we adopt 
a cosmology with 0,m= 0.3, Ha = 0.7 and Ho = 70 km s _1 
Mpc -1 , so that 1" corresponds to 4.22 kpc at the redshift 
of A689. We define spectral index, a, in the sense S tx v~ a . 



mean rate=1 .6821 s 1 



2 OBSERVATIONS AND DATA REDUCTION 

The Chandra observation of A689 (ObsID 10415) was car- 
ried out January 01, 2009. A summary of the cluster's prop- 
erties is presented in Table [T] The observation was taken in 
VFAINT mode, and the source was observed in an ACIS-I 
configuration at the aim point of the 13 chip, with the ACIS 
S2 chip also turned on. 

For the imaging analysis of the cluster we used the 
ciacQ 4.2 software package with cald^I version 4.3.0 and 
followed standard reduction methods. Since our observation 
was telemetered in VFAINT mode additional background 
screening was applied by removing events with significantly 
positive pixels at the border of the 5x5 event islancjf]. We in- 
spected background light curve s of the observatio n following 
the recommendations given in iMarkevitch et all l|2003h . to 
search for possible background fluctuations. The light curve 
was cleaned by 3<r clipping and periods with count rates 
>20% above the mean rate were rejected, curve with re- 
jected periods showed in red. The final level-2 event file had 
a total cleaned exposure time of 13.862 ks. 

As discussed in Sect. 4, there is a bright point source at 
the center of the extended cluster emission which is affected 
by pileup. For the analysis of this source we followed the 
same reduction method, with the exception that VFAINT 
cleaning was not applied. Applying VFAINT cleaning leads 
to incorrect rejection of piled-up events, introducing arti- 
facts in the data. 



1 See http://cxc.harvard.edu/ciao/ 

2 See http:// cxc . harvard . edu/caldb/ 

3 See http://cxc.harvard.edu/ciao/why/aciscleanvf.html 
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Figure 1. Background light curve for the observation of A689 
in the 0.3-12.0 keV band. The CCD on which the cluster falls 
(ACIS-I3) and all point sources are excluded. The red bands show 
periods excluded by the Good Time Interval (GTI) flic. 



source 
background 



J2 le-07 



£1 

8 1e-08 



— f 1 — * — ^^^* h 

10 100 1000 

radius (arcsec) 

Figure 3. Exposure-corrected radial surface brightness profiles 
(0.7-2.0 keV band) of the cluster (red) and the blank-sky back- 
ground (blue), with the green line representing the radius beyond 
which no significant cluster emission is detected. 



3 X-RAY CLUSTER ANALYSIS 

In this section we determine global cluster properties of 
Abell 689. Figure [2] shows a Gaussian smoothed image of 
the cleaned level 2 events file in the 0.7-2.0 keV band (the 
readout streak is removed using the ciao tool acisread- 
corr), with an inset image of the point source which lies at 
the center of the cluster. The extent of the diffuse cluster 
emission was determined by plotting an exposure-corrected 
radial surface brightness profile (Fig [3}, in the 0.7-2.0 keV 
band, of both the observation and blank-sky background 
to determine where the cluster emission is lost against the 
background. Figure [3] demonstrates that the diffuse cluster 
emission is detectable to a radius r w 570" (~ 2.41 Mpc). 
At large radii (r^700") the curves rise due to vignetting 
corrections (larger at larger radii) applied to all the counts, 
whereas in reality each curve contains a component from 
particles that have not been focused by the telescope. 
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Name RA 


Dec 


Z N H ,Gal 


RASS L x ,0.1-2.4kcV 


Abell 689 08 h 37 m 24 s . 70 


+14° 58' 20". 78 


0.279 3.66X10 20 cm" 2 


30.4 xlO 44 erg s" 1 



Table 1. Basic properties. Columns: (1) = Source name; (2) = Right Ascension at J2000 from Chandra; (3) = Declination at J2000 from 
Chandra; (4) = Redshift; (5) = G alactic column density; (6) = Intrinsic X-ray luminosity in the 0.1—2.4 keV band based upon ROSAT 
observations labeling et al.lll99g|) . 
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Figure 2. 0.7-2.0 keV image of A689, smoothed by a Gaussian (cr=1.5 pixels, where 1 pixel = 3.94"), cleaned in VFAINT mode and 
with the readout streak removed. Inset is a zoomed-in unbinned imago of the central point source within A689, cleaned in FAINT mode 
(see § [2j . The inner black circle (r = 208") is the region in which we extract the spectra for our analysis of the cluster emission (see 
§ 13. 2D . the outer black circle (r = 570") represents the detected cluster radius (see §0, and the region between this and the black box 
was used for the local background (see § 13.211 . The white box displays the size of the inset, and the inner white circle (r = 26") shows 
the region excluded due to the central point source. Many point sources are seen in the observation and are excluded from our analysis. 
An extended source on the NE chip can be seen, that is unrelated to A689 and is also excluded from our analysis. 



3.1 Background Subtraction 

In order to take the background of the observation into 
account, appropriate period E Chandra blank-sky back- 
grounds were obtained, processed identically to the cluster, 
and reprojected onto the sky to match the cluster observa- 
tion. We then followed a method similar to that outlined in 
IVikhlinin et al.l (|200fj| ). in order to improve the accuracy of 
the background by applying small adjustments to the base- 
line model. Firstly we corrected for the rate of charged par- 
ticle events, which has a secular and short-term variation 
by as much as 30%. We renormalise the background in the 
9.5-12 keV band, where the Chandra effective area is nearly 
zero and the observed flux is due entirely to the particle 



background events. The renormalisation factor was derived 
by taking the ratio of the observed count rate in the source 
and background observations respectively. The normalised 
spectrum of the blank-sky background is shown in Figure [4] 
over-plotted on the local background for comparison. The 
spectra agree well in the 9.5-12.0 keV band, and across the 
whole spectrum with only slight differences. In addition to 
the particle background, the blank-sky and source obser- 
vations contain differing contributions from the soft X-ray 
background, containing a mixture of the Galactic and geo- 
coronal backgrounds, significant at energies ^1 keV. To take 
into account any difference in this background component 
between the blank-sky and source observations, the spectra 
were subtracted and residuals were modeled in the 0.4-lkeV 
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Figure 4. Comparison of the local (black) and blank-sky (red) 
background spectra, normalised to match in the 9.5-12.0 keV 
band. 



band using an APEC thermal plasma model (|Smith et al.l 
2001). This model was included in the spectral fitting for 
the cluster analysis. As can be seen in Figure [4] this compo- 
nent is very weak in the case of A689. 
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Figure 5. Luminosity- Temperat ure relation of a sample of 115 
clusters of iMaughan et al] jjooi ) fblue open circles). The lumi- 
nosities are measured within [0 < r < ljrsoo and the temperatures 
within [0.15 < r < l.Ojrsnn, in order to minimize the effect of cool 
cores on the derived cluster temperature. Our derived tempera- 
tures for A689 are overplotted for comparison (pink square, green 
triangle, red diamond) (see § 13.21 1. 



3.2 X-ray Cluster Properties 

The analysis of the diffuse X-ray emission allows us to deter- 
mine the X-ray environment surrounding the cluster central 
point source. Throughout this process we excluded the cen- 
tral point source (Fig [2] r^26") and the associated readout 
streaks to avoid contaminating the cluster emission. 

To determine cluster properties we extract spectra out 
to a radius chosen so that the cluster has the maximum 
possible signal-to-noise (SNR). The net number of counts, 
corrected for background, is then 414 (with SNR = 15) 
and the extraction annulus is 26"<r<208" (see Fig[2]), cen- 
tered on the cluster (at a, 8 = 08 h 37 m 24 s . 70, 14° 58' 
20". 78). We fitted the extracted spectrum in xspec with 
an absorbed thermal plasma model (WABSxAPEC) and 
subtracted the background described in § 13.11 We obtain a 
temperature of 13.6^5 3 { 2 keV and a bolometric luminosity 
of L bo i=(10.2±2.9) x 10 44 ergs s -1 . The measured temper- 
ature is far above what we would expect from the lumi- 
nosity. Figure [5] shows the luminosit y-temperature relatio n 
for a sample of 115 galaxy clusters (jMaughan et al] 12008). 
along with the luminosity and temperature derived for A689 
from o ur values above (pink square). As the lMaughan et"al] 
(2008) sample of clusters covers a wide redshift range, the 
luminosities of the clusters were corrected for the expected 
self-similar evolution, given by LxxE(z) -1 , where: 

e{z) = n M0 (i + zf + (i - n M0 - n A )(i + zf + n A . (i) 

The same correction was also applied to the A689 data 
for the plot. Our luminosity was derived within the same 
annular region as the cluster temperature and extrapolated 
both inward and outward in radius between (O-l)rsoo (where 
rsoo represents the radius at which the density of the clus- 
ter is 500 times the critical density of the Universe at that 
redshift) using parameters from a /3-profile fit to the surface 
brightness profile. This takes into account our exclusion of 



the region near the central point source within the cluster, 
an d the extrapolat i on to rsno is in order to compare with 
the lMaughan et all (2008) sample. Our r 5 oo value was deter- 
mined from the te mperature and us i ng the relation between 
I'soo and T given in lVikhlinin et al] (|2006l ). 

The high temperature for A689 could be an indication 
that our observation suffers from background flaring which 
would lead to an overestimate of the cluster temperature. 
However, no evidence is found that the spectrum of the 
blank-sky background differs from that of the local back- 
ground (Figure[4]| , or that the background of the observation 
suffers from periods of flaring (Figure [l}. To investigate the 
sensitivity of the temperature to our choice of background, 
we repeated the analysis using a local background region far 
from the cluster emission (see Fig[2}. We obtained a temper- 
ature of lO.Olg 3 ^ 8 keV. This temperature is again anomalous 
given the luminosity-temperature relation (Fig[5j green tri- 
angle. The luminosity is derived using the same method as 
above, only this time using this new value for the tempera- 
ture to derive rsnn. The spectrum with a local background 
subtracted is shown in Figure [6] We see an excess of pho- 
tons in the ~6-9 keV band, which might have an effect on 
the spectral fit. We perform the same fit using a local back- 
ground, however this time fitting in the 0.6-6.0 keV band to 
ignore these excess high energy photons. We obtain a tem- 
perature of 9.4t2;9 keV. 

Figure [5] shows that the temperature is high for the 
luminosity, or the cluster is X-ray under-luminous. Before 
considering a physical interpretation for the apparently high 
cluster temperature, we investigated possible systematic ef- 
fects from the background subtraction. This was done by 
independently modeling the background. 

We model the background based upon a physical rep- 
resentation of its components. Our model consists of a ther- 
mal plasma APEC model, two power-law components and 
five Gaussian components. The APEC model and one of the 
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Figure 6. Spectrum of the cluster with the local background 
subtracted fitted with an absorbed thermal plasma model, with 
a reduced statistic of x^=l-45 (^=79). 
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Figure 7. Spectrum of the blank-sky background in the source 
region and corresponding the fit (see Sect. 3.2), x2=l-17 (f=585). 



power-law components are convolved with files describing 
the telescope and instrument response (the Auxiliary Re- 
sponse File (ARF) and Redistribution Matrix File (RMF)), 
and are intended to model soft X-ray thermal emission from 
the Galaxy and unresolved X-ray background. The second 
power-law component is not convolved with the ARF as it is 
used to describe the high energy particle background, which 
does not vary with effective area. The five Gaussian com- 
ponents are similarly not convolved, as these are used to 
model line features in the spectrum caused by the fluores- 
cence of material in the telescope and focal plane caused by 
high energy particles. 

As described below and summarized in Tabled the pa- 
rameters of this background model were fit to the blank- 
sky background, local background or high energy cluster 
regions or taken from the literature, in order to build the 
most reliable model possible. We start by modeling blank- 
sky background in the region we used to extract our cluster 
spectrum, with the model outlined above. This allows us to 
place reasonable constraints on the line energy and widths 
of each Gaussian component, and the slope of the uncon- 
volved power-law. We find line features at energies of 1.48, 
1.75, 2.16, 7.48 and 8.29 keV. We also model the convolved 
power-law component, fixing the slope at a value of 1.48 
taken from lHickox fc Markevitchl (|2006T l. The normalisation 
of this power-law component will then be used throughout 
our modeling process, scaled by area where necessary. A 
spectrum of the blank-sky background and a fit using our 
model are shown in Figure [7] 

Next we further constrained model parameters by fit- 
ting our high-energy Gaussian and unconvolved power-law 
components in the 5.0-9.0 keV band within the cluster re- 
gion. At these energies, the emission is dominated by parti- 
cle background. We fix the line energies and widths to those 
found in the blank-sky background and fit for the normal- 
isations. This fit finds a slope of the power-law consistent 
with that found in the blank-sky background. We therefore 
fix the slope of the power-law at r=0.0061, as found for 
the blank-sky background, and fit for the normalisation. We 
finally fit the low energy Gaussians and APEC normalisa- 
tion model parameters in a local background region far from 
the cluster emission. The high energy Gaussians and uncon- 



volved power-law components are frozen at the values found 
in the blank-sky and cluster regions, with the normalisa- 
tions scaled by area. The convolved power-law component 
is frozen at the values found in the blank-sky background 
and the normalisation is scaled by area. The Gaussian fea- 
tures at 1.48, 1.75 and 2.16 keV are frozen at the energies 
and widths found in the blank-sky background. The tem- 
pera ture of the APEC model is frozen at 0.177 keV (taken 
from iHickox fc Markevitchl [20061 ). We note that our APEC 
temperature is not well constrained by our data, but this is 
a weak component. A spectrum of the local background and 
the corresponding fit with the model are shown in Figure [8] 
We now model the cluster with an absorbed thermal 
plasma (WABSxAPEC) model, including the background 
model outlined above. The normalisations of the background 
APEC component and the Gaussians at energies 1.48, 1.75 
and 2.16 keV were scaled by the ratio of the areas from the 
local background region to the cluster region. The normali- 
sations of the fluorescent lines also vary with detector posi- 
tion. To account for this effect in the low energy (<3 keV) 
lines (where we must fit the normalisations in the off-axis 
local background region), we measure the relative change in 
the normalisations of each line between the local background 
and source regions in the blank-sky background data. The 
normalisation of each low energy Gaussian component in the 
fit to the cluster data is scaled for the different detector of 
the cluster region by the relative change in normailsiation 
determined above (in addition to the geometrical scaling for 
size of the extraction region). The unconvolved power-law 
and Gaussians at 7.48 and 8.29 keV are all frozen at the 
values found in the 5.0-9.0 keV cluster region fit. The nor- 
malisation of the convolved power-law is frozen at the value 
found in the blank-sky background. All parameters of the 
model used to describe the background are frozen in the 
corresponding cluster fit, we also freeze the redshift at 0.279 
and the abundance at 0.3. Our fit yields a temperature of 
5.11^'g keV (x^=l-15 (^=79)). We measure a bolometric 
luminosity of Lboi = 1.7xl0 44 erg s _1 . The result is shown 
in Figure [5] (red diamond). The spectrum with the corre- 
sponding fit to the cluster including the background model 
is shown in Figure [9] 
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Component 


Represents 


Parameter 


Value 


Where measured 


Convolved power-law 


Unresolved X-ray bg 


slope 


1.48 


Hickox & Markcvitch (2006) 


normalisation 


4.11xl0 -b 


blank-sky bg 


Unconvolved power-law 


Particle bg 


slope 


0.061 


blank-sky bg 


normalisation 


0.015 


5.0-9.0 keV cluster region 






energy 


1.48 keV 


blank-sky bg 


Gaussian 1 


Al Ka fluorescence 


width 


0.022 keV 


blank- sky bg 






normalisation 


1.82X 10~ 4 


local bg 






energy 


1.75 keV 


blank-sky bg 


Gaussian 2 


Si Ka fluorescence 


width 


0.95 keV 


blank-sky bg 






normalisation 


1.45x 10 -2 


local bg 






energy 


2.16 keV 


blank-sky bg 


Gaussian 3 


Au Map fluorescence 


width 


0.045 keV 


blank-sky bg 






normalisation 


2.24xlO~ J 


local bg 






energy 


7.48 keV 


blank-sky bg 


Gaussian 4 


Ni Ka fluorescence 


width 


0.022 keV 


blank-sky bg 






normalisation 


5.73xlO~ J 


5.0-9.0 keV cluster region 






energy 


8.29 keV 


blank- sky bg 


Gaussian 5 


Cu + Ni fluorescence 


width 


0.168 keV 


blank- sky bg 






normalisation 


4.56x10 


5.0-9.0 keV cluster region 






kT 


0.177 keV 


Hickox & Markcvitch (2006) 


APEC 


Galactic foreground emission 


abundance 


1.0 


solar abundance 


redshift 





Galactic 






normalisation 


2.9xl0~ a 


local bg 



Table 2. Table of the individual model components used to represent the background, with a brief interpretation of each compo- 
nent, individual component parameters, parameter values and where each value is calculated. All normalisations are measured in 
photons/keV/cm 2 /s at 1 keV, and scaled to the cluster region. Blank-sky background parameters were derived in the same region 
and the local background were measured in an area 2.45 times that of the cluster and therefore the normalisations reduced by this factor. 
The low energy Gaussians were also corrected due to the normalisation dependence with position on the detector (see text). 
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Figure 8. Spectrum of the local background and corresponding 
fit (see Sect 3.2). x^=0.79 (i/=120). 



Figure 9. Spectrum of the cluster plus background fit with an 
absorbed thermal plasma model including a background model 
(see Sect 3.2). xt=1.15 (u=79). 



4 THE CENTRAL POINT SOURCE 

The point source is displayed in Figure I10I The presence 
of strong readout streaks indicate that the point source is 
likely to be affected by pile-up. The readout streaks occur as 
X-rays from the source are received during the ACIS paral- 
lel frame transfer, which provides 40^ts exposure per frame 
in each pixel along the streak. We detail our analysis of the 
point source and estimate of the pileup fraction in the fol- 
lowing section. 



4.1 X-ray Analysis of the Point Source 

As a first test of the predicted pile-up, we compared the im- 
age of the point source to the Chandra Point Spread Func- 
tion (PSF). We made use of ciao tool mkpsf to create an 
image of the o n-axis PSF following the method outlined in 
iDonato et all l|2003r i. This consists of merging 7 different 
monochromatic PSFs chosen and weighted on the basis of 
the source energy spectrum between 0.3 and 8 keV. This 
method can be summarized as follows: 

(i) We first extract the energy distributions of the pho- 
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Figure 10. Chandra image of A689 showing the regions used for 
extracting spectra of the readout streak (inner rectangles) and the 
corresponding background regions for the readout streak (outer 
rectangles). 
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Figure 11. Surface brightness profiles of the point source (red 
squares) and the composite PSF (blue circles) in the 0.3-8.0 keV 
band, normalised to agree in the 2.46-4.92 arcsec radii region. 

tons from a circular region centered on the peak brightness 
of the source with a radii of 2.5". 

(ii) We choose seven discrete energy values at which to 
creating each PSF, with the number of counts at each energy 
corresponding to that PSF's 'weight'. 

(iii) Using mkpsf we create seven monochromatic PSFs at 
the position of the point source on the detector and co-added 
them. Each PSF is weighted by its relative normalisation 
(found in the previous step). 

Once we obtained the composite PSF, we normalised it 
to the counts within an annulus (inner and outer radius 2.46 
and 4.92 arcsec respectively) in order to avoid any effects 
of pile-up in the central region. We then compare surface 
brightness profiles of the point source and PSF to look for 
evidence of pile-up in the core of the point source image. 
Figure [TT] shows the radial surface brightness profiles of the 
point source (red) and the composite PSF (blue). We find 
that the point source and PSF agree well in the wings of the 



PSF (> 2.46") and that there is an excess in the PSF surface 
brightness above that of the point source at the peak of the 
source. This is consistent with the flattening of the source 
profile relative to that of the PSF due to pileup in the core. 
The PSF then gives an estimate of the non piled up count 
rate. Given this count rate we use PIMMS3 to estimate a 
pile-up fraction of 65%. 

We also compute a second estimate of the core count 
rate using the ACIS readout streak. By fitting a model 
to the spectrum extracted from the readout streak we can 
compare this to a spectrum extracted in the core and fit- 
te d using a pileup mode l. We follow the method outlined 
in iMarshall et all |2005) in order to correct the exposure 
time of the readout streak spectrum. For an observation 
of live time ti ive , a section of the readout streak that is 
8 3 arcsec long accumulates an exposure time of t s = 4 
xl0 _5 ti ivc 6> s /(t/6> x ) s, where 6> x = 0.492" is the angular 
size of an ACIS pixel. For our observation, tii vo = 13.862 
ks and the frame time parameter tf = 3.1 s, giving t a 
= 165s in a streak segment that is 454" long. Figure 1101 
shows the regions used for extracting spectra of the readout 
and an adjacent background region. This choice of back- 
ground region ensures the cluster emission is subtracted 
from the readout str eak spectrum. Using the sherpa package 
( Rcfsc tal et al.l [2009) we fitted an absorbed 1-D power-law 
model (WABSx POWER-LAW) to the extracted spectrum 
of the readout streak. We obtain fit parameters for the pho- 
ton index = 2.33±° ?/ and a normalisation of 0.0033±0.0005 
photons keV" 1 cm -2 s" 1 (x 2 =0.34 (z^=63)). We then ex- 
tract a spectrum of the point source in a region of radius 
5", and subtracted the same background as for the readout 
streak. We once again fitted an absorbed power-law model, 
including this time a pileup model (jdpileup). We obtain fit 
parameters for the photon index = 2.22^'^, consistent to 
that found from the readout streak, and a normalisation of 
0.0015±0.0001 photons keV" 1 cm -2 s _1 (xl=l-6 (^=119)). 
The extracted spectrum and corresponding fit is shown in 
Figure The pileup fraction is estimated to be 60%, which 
is consistent with that found using the PSF count rate. The 
normalisation found in the fit can be converted to an X-ray 
flux density for this source, fi kcv=0.99±0.07 fiJy. 

4.2 Optical Observations 

Abell 689 was observed with the Hubble Space Telescope 
(HST) with the F606W filter (V-band) on January 20, 2008. 
Marking the position of the peak X-ray emission of the point 
source on the HST image, we find that this corresponds to 
an object resembling an active nucleus in a relatively bright 
galaxy (Fig [13}. We searched the SDSS DR7 archive for in- 
formation on the spectral properties of this object. At the co- 
ordinates of the X-ray point source (SDSS coordinates of a, 5 
= 08 h 37 m 24.7, 14° 58' 19". 8) we find a blue object with a 
corresponding relatively featureless spectrum (Fig |14jl . From 
SDSS we quote an r-band magnitude of 17.18. The spectrum 
resembles that of a BL Lac object, a type of AGN orientated 
such that the relativistic jet is closely aligned to the line 
of sight. From the H and K lines in the spectrum (dotted 
green lines in Fig |14|l . thought to be from the host galaxy, 

4 http:/ /cxc. harvard.edu/toolkit/pimms.jsp 
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Energy (keV) 



Figure 14. Spectra of the BL Lac from SDSS (SDSS 
J083724. 71+145819. 8). The vertical green lines represent the H 
Figure 12. Spectra of the point source fitted with an absorbed and K lines thought to be from the host galaxy and the green line 

power-law, including a pileup model. y^=\.§ (^=119). at the bottom represents the error spectrum. 
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Figure 13. HST image of the point source with a cross marking 
the position of the peak of the X-ray emission. 

the redshift is determined to be z =0.279, consistent w ith the 
redshift assigned to the cluster l|Collins et al.lll995r i. Using 
the HST observation we measure an optical flux for the BL 
Lac of f °=112mJy. 



4.3 Radio Observations 

Archival radio observations of Abell 689 are available, allow- 
ing us to determine the radio spectral index, a. We obtained 
8.46 GHz data taken in March 1998 from the VLA archive 
which we mapped using AIPS. The source is unresolved at 
8.46 GHz, and we measure a flux density of 18.6T0.27 mjy. 
We also obtained a 1.4 GHz radio image from the FIRST sur- 
vey, from which we measure a flux of 62.7±0.25 mjy. From 
these data we obtain a spectral index of a r ~0.67±0.01. We 
note that a slight angular extension in the FIRST survey 
suggests that the 8.46 GHz image may be missing some flux 
density. 



5 DISCUSSION 

5.1 The ICM properties of Abell 689 

We have derived the ICM properties of A689 using three 
methods of background subtraction. Table|3]shows the spec- 
tral properties of the ICM for each of the background treat- 
ments we employ with our favored method coming from a 
physically motivated model of the background components. 
Through detailed modeling of the local and blank-sky back- 
grounds for A689, and including this background model in 
a spectral fit to the cluster, we have determined a tempera- 
ture and luminosity for A689 of T = 5.1 jl^'g keV and Lboi = 
3.3xl0 44 erg s _1 . Plotting these values on the luminosity- 
temperature plot (Fig O and comp aring to the large X- 
ray sample of iMaughan et all |200ct ). we find that A689 is 
observed to be at the edge of the observed scatter in the 
luminosity-temperature relation. This suggests that either 
the temperature of the ICM has been enhanced or suppres- 
sion of the luminosity has occurred. 

It has been shown that systems that host a radio source 
are likely to have hig her temperatures at a given luminosity. 
ICroston et alj (|2005l ) sho wed that this is the case fo r radio- 
loud galaxy groups, and iMagliocchetti fc Briiggenl (|2007f ) 
showed that for clusters that host a radio source there is a 
departure from the typical luminosity-temper ature relation, 
partic ularly in the case of low mass systems. ICroston et al.l 
(2005) also showed, through analysis of Chandra and XMM- 
Newton observations, evidence for radio-source interaction 
with the surrounding gas for many of the radio-loud groups. 
A similar process could be occurring within A689, which 
has a confirmed radio source at the center of the cluster. We 
note that more detailed X-ray and radio observations would 
to be needed in order to test any interaction between the 
BL Lac and ICM. 

The other possible explanation for the offset of radio- 
loud systems from the luminosity-temperature relation is 
the suppression of the luminosity, which could be caused 
by displacement of large amounts of gas due to the inter- 
action of the radio source with the ICM. The interaction 
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10.0+«3 8 

cc 1+2.2 
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6.2+1.4 
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1.19 
1.45 
1.15 


74 
79 
79 



Table 3. Table listing the derived spectral properties of A689 for each of the background treatments we employ in our analysis, t 
indicates our favored method for determining the cluster properties of A689. 



of the radio source with the ICM will cause an increase in 
the entropy of the local ICM. This higher entropy gas will 
be displaced so that it is in entropy equilibrium w i th the 
surrounding gas. However, iMagliocchetti &: Brugger] (|2007l ) 
showed that there is a correlation between the radio luminos- 
ity and the heat input required to produce the observed tem- 
perature increment in clusters hosting radio sources. They 
note that this correlation favors an enhanced temperature 
scenario caused by the radio galaxy induced heating. 

5.2 Comparison with the BCS 

A689 was noted in the BCS as having a significant fraction 
of its flux coming from embedded point sources. We have 
confirmed that A689 contains a point source at the center 
of the cluster and is that of a BL Lac object. We stated 
that the measured X-ray luminosity for A689 as quoted in 
the BCS (Lo.i_2.4kov = 3xl0 45 erg s" 1 , see § Q}, is the 
third brightest in the BCS. From our follow up observa- 
tion with Chandra we calculate the luminosity and compare 
with the BCS value. The luminosities in the BCS are cal- 
culated within a standard radius of 1.43 Mpc, which at the 
redshift of A689 corresponds to a radius of 338". We there- 
fore employ the same method as in Sect 13.21 and integrate 
under a beta model fitted to the derived surface brightness 
profile and extrapolate inward and outward from 26-206" to 
0-338". The unabsorbed flux in the 0.1 - 2.4 keV band (ob- 
served frame) was fo.i-2.4,keV — 5.8xl0 -13 erg s _1 cm -2 . 
After k-correction the X-ray luminosity in the 0.1 - 2.4 keV 
band (rest frame) was Lo.i_2.4,k C v = 2.8xl0 44 erg s _1 . Note 
that we assume an Ho of 50 for consistency with the BCS 
catalog. This value is ~10 times lower than that quoted in 
the BCS, and A689 is now ranked 110th out of 201 in lumi- 
nosity. 



5.3 Classifying the BL Lac 

BL Lac objects may be split into 'High-energy peak BL Lacs' 
(HBL) and 'Low-energy peak BL Lacs' (LBL), for objects 
which emit most of their synchrotron power at high (UV- 
soft-X-ray) or low (far-IR, n ear-IR) frequencies respectively 
|Padovani fc Giommil 1 19951 ). HBL and LBL objects have 
radio-to-X-ray spectral indies of a rx ^0.75 and a ra ^0.75, 
respectively. We calculate a radio-to-X-ray spectral index 
for the BL Lac in A689 of a ra =0-58±0.04. From this value 
we classify our BL Lac as an HBL type. We al so compared 
our BL Lac with those of iFossati et al.l (| 19981 ). who inves- 
tigated the properties of large samples of BL lacs at ra- 
dio to 7-ray wavelengths. Our value of a rx —0.58 falls into 
the region 0.35^a r a:^0.7, dominated by X-ray selected BL 



Lacs (XBL). This is consistent with the X-ray selection of 
this cluster. Using the measured HST flux, we also calculate 
a ro =0.50 and aw =0.77. T hese values are not atypical for 
BL Lac objects (F i gure 7 in lWorrall et al.ll 19991 ). 

iDonato et al.l (|2003i) tried to determine whether HBLs 
and LBLs were characterized by different environments. 
They found that of 5 sources exhibiting diffuse X-ray emis- 
sion that 4 were HBLs and 1 was an LBL. The BL Lac 
in A689 continues this trend, as it appears to be an HBL 
embedded within a cluster environment. 

5.4 Evidence for Inverse-Compton Emission 

Our initial analysis of the extended emission in A689 yielded 
temperature estimates that were significantly higher than 
expected based upon its X-ray luminosity. We also found 
evidence for an hard excess of X-ray photons in the 6.0-9.0 
keV band of the cluster spectrum. Before assigning a physi- 
cal cause we must be careful to eliminate systematic effects 
in the background subtraction, as underestimating the par- 
ticle background will give a hard excess. This is unlikely to 
be the case here as a high temperature was obtained when 
either a blank-sky and local background was used. 

When using our background model instead of subtract- 
ing a background spectrum, the cluster temperature was 
more consistent with its luminosity (Fig [5]) . This appears 
to be due to the unconvolved power law fitting a hard ex- 
cess in the cluster region. In order to assess what impact 
the high-energy particle component has on the cluster tem- 
perature, we varied the normalisation of the unconvolved 
power-law component of our background model (sect. 3.2) 
within its la errors. The temperature ranged from 3.86 to 
8.51 keV. Thus small changes in the particle background 
have a significant influence on the cluster temperature. Our 
power-law component was derived within the cluster region 
in the 5.0-9.0 keV band, so our background model may be 
removing a physical hard X-ray component associated with 
the cluster. We assess the possible properties of such a com- 
ponent by re-deriving the normalisation of the unconvolved 
power-law component in the local background region, scaling 
this value by the ratio of the areas, and using this value for 
the unconvolved power-law normalisation in the background 
model for our cluster fit. We find a power-law normalisation 
of 0.0129 photons keV" 1 cm" 2 s" 1 (as opposed to 0.0149 
photons keV" 1 cm -2 s _1 in the cluster region). Using this 
value and re-fitting (Fig|15[). we find a cluster temperature 
of 14.3+g 5 4 6 keV. This suggests that the hard component is 
spat ially associated w i th the cluster. 

iBonamente et ail |2007l ) find a similar excess of X-ray 
emission in the cluster A3112, which is known to have a 
central AGN. It is argued that this excess may be due to 
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Figure 15. Spectrum of the cluster including our physically mo- 
tivated model, with the normalisation of the unconvolved power- 
law component derived in the local background region. x^=\.Y! 
(v=79) 



emission of a non-thermal component. Relativistic electrons 
in the intergalactic medium will cause CMB photons to scat- 
ter into the X-ray band (inverse Compton scattering). The 
same process could occur in A689, with the relativistic parti- 
cles responsible for the inverse Compton scattering provided 
by the jets of the BL Lac. 

In order to test the assumption of inverse-Compton 
emission, we add a convolved power-law component with 
a— 1.5 (appropriate for modeling IC emission from aged elec- 
trons). We fit for the normalisation of this added power- 
law component and freeze all other parameters of the back- 
ground and cluster model. We note that we use the normali- 
sation of the unconvolved power-law found in the local back- 
ground region as found above (a value of 0.0129). When fit, 
the x 2 increases slightly but the fit is still acceptable at the 
95% confidence level. From the additional power-law com- 
ponent we measure a 1 keV flux density of ~7 njy. If the 
X-ray emission is from scattering of the CMB by an aged 
population of electrons of power-law number index 4.0, we 
can determine what the implied synchrotron emission would 
be at 1.4 GHz, given plausible magnetic fiel ds. Clusters typ- 
ically have magnetic fields of a few fiG (jCarilli fc Taylor! 
120021) . The NVSS survey would have detected a flux den- 
sity of ~10 mjy over the entire cluster. Adopting a 1.4 GHz 
flux density limit of 10 mjy, we require a magnetic field of 
B ^ 2 /j,G over the cluster to avoid over-predicting radio 
emission. We also calculate a value for the minimum-energy 
magnetic field, BE_min, that would give a 1.4 GHz flux den- 
sity of Si. 4 GHz = 10 mjy, which in this case is BE_min = 7.5 
fxG. For a flux density of Si. 4 ghz = 5 mjy, we would re- 
quire a magnetic field of 1.6 /iG and BE_min = 6.5 fiG. Note 
that B oc S1.4 ghz 1/1+ " and B E _ mi n <x S1.4 ghz 1/3+q . We 
conclude that the excess X-ray emission can be attributed 
to inverse-Compton scattering without over-predicting the 
radio emission if the magnetic field strength is in a range 
typical of clusters and within a factor of a few of the mini- 
mum energy value. 



6 CONCLUSIONS 

We have used a 14 ks Chandra observation of the galaxy 
cluster A689 in order to determine the nature of the clus- 
ter's point source contamination and to analyze the cluster 
properties excluding the central point source. Our main con- 
clusions are as follows. 

1. Using background subtraction of both local and blank- 
sky backgrounds, we obtain temperatures which are high 
relative to the luminosity-temperature relation. 

2. We construct a physically motivated model for the back- 
ground and include this model in a fit to the cluster spec- 
trum. If the particle background in the cluster is allowed 
to exceed that in the local and blank-sky backgrounds we 
obtain a temperature of 5.lj^'g keV. However, there is no 
reason for there to be a higher particle rate in the specific 
region of the CCD in which the cluster lies. A hard excess 
needed to bring the temperature to a reasonable value must 
have a different origin. 

3. We confirm the presence of a point source within A689 
as suspected in the BCS. When excluding the point source 
and using our derive background model we find a luminosity 
of Ln.i_2.4,k c v = 2.8xl0 44 erg s" 1 , a value ~10 times lower 
than quoted in the BCS. 

4. From the X-ray analysis of the point source we find a 
"flat-topped" point source with a pileup fraction of ~ 60%. 

5. Optical observations of the cluster from SDSS and HST 
lead us to conclude that the point source is a BL Lac type 
AGN. 

6. We classify the BL Lac as an 'High-energy peak BL Lac' 
with « ra; =0.58±0.04. 

7. We interpret the hard X-ray excess needed to bring 
the cluster temperature to a reasonable value as inverse- 
Compton emission from aged electrons that may have been 
transported into the cluster from the BL Lac. 

We have shown here not only the importance of re- 
solving and excluding point sources in cluster observations, 
but also the effect these point sources can have when de- 
termining the ICM properties of galaxy clusters. The de- 
tailed analysis we have performed here may not however be 
suitable for all clusters as it is unclear whether this anal- 
ysis can be performed at higher redshifts. Separating the 
point source and cluster emissions becomes increasingly dif- 
ficult at high redshifts, however Chandra has proved capa- 
ble at resolving point so urce emission in clusters to z~2 
fe.g. lAndreon et al.ll2009h . Resolving point source and clus- 
ter emission out to high redshift becomes more important 
with redshift due to the evolut ion of the number d ensity 
of point sources within clusters (|Galametz et al.ll2009l ). The 
area used to model the background components associated 
with the high energy cluster regions will decrease spatially 
with redshift, and separating the between thermal and in- 
verse Compton emi ssions at higher red shift will become in- 
creasing ly difficult (|Fabian et al.ll2003l ). due to the increase 
in the energy density of the CMB with redshift. 
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